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Based on the analysis of the measurement data of angle-resolved photoemission spectroscopy 
(ARPES) and optics, we show that the charge transfer gap is significantly smaller than the optical 
one and is reduced by doping in electron doped cuprate superconductors. This leads to a strong 
charge fluctuation between the Zhang-Rice singlet and the upper Hubbard bands. The basic model 
for describing this system is a hybridized two-band t-J model. In the symmetric limit where the 
corresponding intra- and inter-band hopping integrals are equal to each other, this two-band model 
is equivalent to the Hubbard model with an antiferromagnetic exchange interaction (i.e. the t-U-J 
model). The mean-field result of the t-U-J model gives a good account for the doping evolution of 
the Fermi surface and the staggered magnetization. 



The evolution of the Fermi surface and the Mott insu- 
lating gap with hole or electron doping is a central issue in 
elucidating the mechanism of high-Tc superconductivity. 
In the hole doped case, a small Fermi surface arc appears 
first near (7r/2, 7r/2) and then extends towards (tt, 0) and 
(0, tt) with increasing doping. In contrast, in the electron 
dope case, electrons are first doped into the upper Hub- 
bard band (Cu 3d^° band) near (tt ,0) and equivalent 
points. With further doping but still in the antiferro- 
magnetic phase, in-gap spectral weight develops below 
the Fermi level. These in-gap states move upwards and 
eventually form a hole-like Fermi surface pocket around 
(7r/2 , 7r/2) 1]. In the heavily overdoped sample, these 
two Fermi pockets merge together and form a large Fermi 
surface with a volume satisfying the Luttinger theorem. 

The peculiar doping dependence of the Fermi surface 
topology in electron-doped cuprates is a manifestation of 
correlation effects. To understand the physics behind, 
much of the theoretical study has been carried out with 
the one-band Hubbard model [3, y, lj|- ^^ ^^^^ model, a 
metallic band is split into two effective bands, namely 
the upper and lower Hubbard bands, by a correlation 
energy U that represents the energy cost for a site to 
be doubly occupied. Under the mean-field approxima- 
tion, this model gives a good account for the experimen- 
tal data if U is assumed to fall strongly with doping. 
However, this strong reduction of U by doping is not 
usually expected [5|. The nominal Hubbard [/-term could 
arise either from the on-site Coulomb repulsion between 
two electrons in a Cu 3d.j.2_y2 orbital or from the charge 
transfer (CT) gap between O 2p and Cu Sd^" bands. The 
on-site Coulomb repulsion in a Cu 3dj.2_y2 states in high- 
Tc cuprates is generally larger than 5 eV. The CT gap 
in electron doped cuprates quoted from the optical data 
is also quite large (~ 1.5 eV)[a, 0|- It seems that in nei- 
ther case U can be dramatically suppressed by only 15% 
doping. 



An alternative interpretation to the two Fermi pockets 
is based on the notion of band folding induced by the an- 
tiferromagnetic interaction [3, Q. This interpretation is 
consistent with the measurement data in the overdoped 
regime {x > 0.15). However, in the low doping anti- 
ferromagnetic phase, it breaks down. The band folding 
assumes implicitly a band with large Fermi surface exists 
and it is the antiferromagnetic interactions between the 
hot spots split this band into a conduction electron and 
a shadow hole band. However, in the antiferromagnetic 
phase at low doping, these bands with the folding gap 
at the hot spots were not observed and the band near 
(7r/2,7r/2) is well below (7r,0).[l[ Furthermore, the anti- 
ferromagnetic interaction is too small to account for the 
energy splitting between the lower and upper CT bands 
at least in the low doping limit. 

To resolve the above problems, it is important to un- 
derstand correctly how the hole-like Fermi pockets de- 
velop with doping. In a nominally undoped Nd2Cu04, 
a dispersive band structure is observed by ARPES at 
roughly 1.2 eV below the chemical potential. As shown 
in Ref. [1] , the energy-momentum dispersion of this spec- 
tral peak behaves almost the same as the lower CT band 
observed in Ca2Cu02Cl2, except in the latter case the 
band lies at only ~ 0.7 eV below the chemical potential. 
This suggests that these two bands have the same physi- 
cal origin. The difference is probably due to the intrinsic 
doping and the Fermi energy is pinned near the bottom 
of conduction band (i.e. Cu 3d^" band) in Nd2Cu04 and 
near the top of the valence band {i.e. Zhang- Rice singlet 
bandflU) in Ca2Cu02Cl2. 

Doping electrons into Nd2Cu04 results in a spectral 
weight transfer from the main spectral peak at ~ 1.2eV 
to a "in-gap" state. This in-gap state first appears as a 
week low energy "foot" at ^ 0.5 eV below the Fermi level 
£p along the zone diagonal in the undoped Nd2Cu04 
(Fig. [T|) . It moves towards the Fermi level with doping 
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Figure 1: ARPES spectra near (a) the nodal and (b) antinodal 
regions, reproduced from the data pubhshed in Ref. ]2^. (c) 
The infrared conductivity reproduced from the data pubhshed 
in Ref. [Tj. The insets of (b) and (c) illustrate the indirect 
and direct CT gaps. 



and becomes a broad hump just below the Fermi level at 
optimal doping. The hole Fermi pockets observed at high 
doping originate from these in-gap states. In contrast, 
the states near (tt, 0) resides at ep as they are derived 
from the bottom of the upper Hubbard band. The fact 
that the broad maximum is slightly below ep is caused 
by the Frank-Condon broadening as discussed below. 

It should be pointed out that, same as for the dis- 
persive high energy band, the in-gap states behave simi- 
larly as the low energy coherent states observed in hole- 
doped Ca2Cu02Cl2[9[. Near half filling, the in-gap state 
in Nd2Cu04 lies also at ~ 0.7 cV above the high en- 
ergy spectral peak. This suggests that, similar as in hole 
doped materials, the high energy hump structure in the 
spectra results from the Franck-Condon broadening and 
the in-gap states are the true quasiparticle excitations lo- 
cated at the top of the lower CT bandQ. At half filling, 
the in-gap state is not observed because its quasiparticle 
weight is vanishingly small [10|. 

The spectral weight transfer induced by doping has 
also been observed in the optical measurements (Fig.[T]). 
At zero doping, the optical CT gap appears at ~ 1.5 
eV. Upon doping, a mid-infrared conductivity peak de- 
velops. This mid-infrared peak appears at ^ 0.5 eV at 
low doping [y, |7|, and then moves towards zero energy 
with increasing doping. The doping dependence of the 
mid-infrared peak is consistent with the doping evolution 
of the "in-gap" states observed by ARPES. It suggests 
that the mid-infrared peak results mainly from the opti- 
cal transition between the "in-gap" states and the upper 
Hubbard band. The polaron effect may ako have some 
contribution to these mid-infrared peaks. ,1^ 

The above discussion indicates that the true CT gap, 
measured as the minimum excitation energy between the 
lower and upper Hubbard bands, is only 0.5 eV at half 
filling, much lower than the optically measured CT gap, 
which is usually believed to be about 1.5 eV. This dif- 



ference between the true quasiparticle gap that deter- 
mines the transport and thermodynamics and optically 
measured CT gaps has also been found in hole doped 
materials. For La2Cu04, Ono et aZ. 13] found recently 
that the CT gap obtained from the high temperature be- 
havior of the Hall coefficient is only 0.89 eV, while the 
corresponding optical CT gap is about 2 eV. This means 
that the optical CT gap, which is generally determined 
from the peak energy of the optical absorption, does not 
correspond to the true gap between the two bands in 
high-Tc oxides. The indirect nature of the gap (insets 
of Fig. [1]) and the Frank-Condon effect lead to the over- 
estimate of the gap by optics. It also means that the 
charge fluctuation in high-Tc materials is much stronger 
than usually believed and should be fully considered in 
the construction of the basic model of high-Tc supercon- 
ductivity [ii,[il[il. 

The doping dependence of low energy peaks observed 
by both APRES and optics indicates that there is a gap 
closing with doping. This gap closing may result from the 
Coulomb repulsion between O 2p and Cu 3d electrons. 
Doping electrons will increase the occupation number of 
Cu 3d states, which in turn will add an effective potential 
to the O 2p states and raise their energy level. If Upd is 
the energy of the Coulomb interaction between neighbor- 
ing O and Cu ions, then the change in the O 2p energy 
level will be Sep « +2xUpd, where x is the doping con- 
centration and the factor 2 appears since each O has two 
Cu neighbors. Upd is generally estimated to be of order 
1-2 eV. Thus a 15% doping of electrons would reduce the 
CT gap by 0.3-0.6 eV, within the range of experimen- 
tally observed gap reduction. Furthermore, the electro- 
static screening induced by doping can reduce the on-site 
Coulomb interaction of Cu 3d^2_y2 electrons. This can 
also reduce the gap between the O 2p and the upper 
Hubbard bands. 
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Figure 2: Fermi surface density map at different dopings x, 
obtained by integrating the spectral function from -40 to 20 
meV around the Fermi level, for the t-U-J model. 



Now let us consider how to characterize the low energy 
charge and spin dynamics of the system. For simplicity, 
we focus on the electronic structure and leave the addi- 
tional clcctron-phonon interaction effect for future study. 
The lower CT band behaves similarly as the Zhang-Rice 
singlet band. Thus if its charge fluctuation with the up- 
per Hubbard band is ignored, this band should be de- 
scribed by an effective one-band t-J model[ll|. Simi- 
larly, the upper Hubbard band should also be described 
by an effective one-band t-J model if there is no charge 
fluctuation. However, in the case the hybridization or 
charge transfer between these bands is important, it can 
be shown from a three-band model that these two t-J 
models should be combined together and replaced by the 
following hybridized two-band t-J model[l7| 
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where hi, Ci and dia are the annihilation operators of a 
Zhang- Rice singlet hole, a doubly occupied dx2_y2 state 
(doublon), and a pure Cu^+ spin, respectively. At each 
site, these three states cannot coexist and the correspond- 
ing number operators should satisfy the constraint 
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The difference between the number of doubly occupied 
(ia;2_j,2 states and Zhang- Rice singlet holes is the doping 

concentration of electrons, (ejei — h\hi) = x. 

In Eq. (HI), Sj = d\adi/2 is the spin operator and a is 
the Pauli matrix, e^ and £h are the excitation energies of 
a doublon and a Zhang- Rice singlet, respectively, t^j and 
t^ are the hopping integrals of the upper Hubbard and 
Zhang- Rice singlet bands. In Eq. ([1]), if e/i :2> £e > 0, 
then {h\hi) « and H simply reduces to the one-band 
t-J model of doubly occupied electrons in the doublon- 
spinon representation. On the other hand, if Eg ^ e^ > 
, then {e\ei) ~ and H becomes simply the one-band t-J 
model of Zhang-Rice singlets in the holon-spinon repre- 
sentation. The tij term describes the hybridization be- 
tween the upper Hubbard and Zhang-Rice singlets. The 
last term results from the Coulomb repulsion between a 
Cu 3dj;2_j,2 and its neighboring O 2px^y electrons. Vpd is 
proportional to the Coulomb repulsion between Cu and 
O ions Upd- 

The above Hamiltonian can be simplified if i? = t^- — 
tij. In this case, by using the holon-doublon represen- 
tation of an electron operator Cia — ah\dicr + eid]-^, 
and taking a mean-field approximation for the Vpd-term, 
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+ Sh- AVpdiielei) + 
(hlhi)). In electron doped materials, as the induced hole 



concentration is very small, (hlhi) <C (e^Ci) « x, we have 
U ~ ee + £h — 4:xVpd- It should be emphasized that the 
spin exchange term in Eq. ([3]) is not a derivative of the 
one-band Hubbard model in the strong coupling limit. It 
actually arises from the antiferromagnetic superexchange 
interaction between two undoped Cu^"*" spins via an O 2p 
orbital. This term, as shown in Ref. [l8|, can enhance 
strongly the superconducting pairing potential. 

The t-U-J model defined by Eq. ([3|) is obtained by 
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t" 



tij . This is a strong approximation 
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which may not be fully satisfied in real materials. Never- 
theless, we believe that this simplified model still catches 
qualitatively the low energy physics of high-Tc cuprates. 
It has already been used, as an extension of either the 
Hubbard or the t-J model, to explore physical proper- 
ties of strongly correlated systems, such as the gossamer 
superconductivity. [19] 

The above Hamiltonian reveals two features about the 
effective Hubbard interaction. First, U is determined by 
the CT gap[14i|, not the Coulomb interaction between two 
electrons in a Cu ?id^-2_y2 orbital. It is in the intermediate 
or even weak coupling regime, rather than the strong 
coupling limit as usually believed. Second, U is doping 
dependent. It drops with doping. These are in fact the 
two key features that are needed in order to explain the 
experimental results with the Hubbard model @, S Ij] ■ 

We have calculated the single-particle spectral function 
and the staggered magnetization for the t-U-J model us- 
ing the mean-field approximation. In the calculation, tij 
are parameterized by the first, second and third near- 
est neighbor hopping integrals (f , i', t"). The parameters 
used are t = 0.326 eV, t' = -0.25i, t" = 0.15i, J = 0.3i 
eV, ee + £h = 4t and Vpd — 2.7i. 

Fig. [5] shows the intensity plot of the spectral function 
at the Fermi level. The doping evolution of the Fermi sur- 
face agrees with the ARPES measurements. [l| It is also 
consistent with the mean-field calculation of the Hubbard 
model by Kusko et al., while our calculation has the same 
shortcoming of the mean field calculation in providing too 
large band width. [2;] The difference is that, in our calcu- 
lation, the Hubbard interaction U is not an adjustable 
parameter of doping. It decreases almost linearly with 
doping. For the parameters given above, U ~ 4t— lO.Sxi. 
Whereas in the calculation of Kusko et al. [2| , U is deter- 
mined by assuming the mean-field energy gap to be equal 
to the experimentally observed value of the "pseudogap" . 

Fig. [3] shows the theoretical result of the stag- 
gered magnetization m — {{ui^ — nii))/2. The sim- 
ple mean-field result agrees well with the experimental 
data[20, [2l|, [22, [23|, especially in the low doping range. 
It is also consistent qualitatively with other theoretical 
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Figure 3: Comparison of the mean field result (open circles) 
of the staggered magnetization m as a function of doping with 
the experimental data (solid circles) 23] . The theoretical data 
obtained by Yan et al. [2J| and by Yuan et al. [3 are also shown 
for comparison. 



calculations p, [23|- 'ti decreases almost linearly at low 
doping. However, it shows a fast drop above ^ 0.14, 
when the lower Zhang-Rice singlet holes begin to emerge 
above the Fermi surface. This abrupt change of m is an 
indication of a significant renornialization of the Fermi 



surface. It may result from the quantum critical fluctua- 
tion as suggested in Ref. [2^. m does not vanish above 
the optimal doping, this is probably due to the mean-field 
approximation. 

In conclusion, we have shown that the minimal CT 
gap is much smaller than the optical gap and the charge 
fluctuation between the Zhang-Rice singlet and the up- 
per Hubbard bands is strong in electron doped copper 
oxides. The low-lying excitations of the system are gov- 
erned by the hybridized two-band t-J model defined by 
Eq. ([T]) or approximately by the t-U-J model defined by 
Eq. ([3]) . This conclusion is drawn based on the analysis of 
electron doped materials. However, we believe it can be 
also applied to hole doped cuprate superconductors. Our 
mean-field calculation for the t-U-J model gives a good 
account for the doping evolution of the Fermi surface as 
well as the staggered magnetization. It sheds light on the 
further understanding of high-Tc superconductivity. 
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